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1 Introduction 



There is a standard framework for computing Green's functions in supersymmetric J\f = 4 

Yang-Mills theory using AdS/CFT duality [H 13 G9 01 El d . In ®' thiS WaS carried out 
explicitly over the entire to, k plane for an operator coupling to a scalar field. The specific goal 
was to compute the real static (to = 0) Green's function both directly and by integrating over 
to. The result was then Fourier transformed to a spatial Green's function and compared with 
a computation in next-to-leading order QCD perturbation theory [10] , see also [TT| IT2] . 

The previous calculation was carried out in a scalefree theory. The purpose of this article 
is to see what effects arise from Aqcd> i-e., scale invariance breaking in SU(iV c ) Yang-Mills 
theory. Physically this manifests itself in the existence of glueballs and a first order phase 
transition^ The predictions are quantitative and it is suggested that their comparison with 
lattice Monte Carlo data may give a better handle on verifying the strongly interacting 
nature of QCD matter in the temperature range T C ...WT C . There is definite phenomenological 
evidence for this (magnitude of pressure, small shear viscosity, jet quenching), but no definite 
nonperturbative test [13] • In fact, how can the system be strongly interacting while so many 
lattice results [15] agree with ideal gas? One concrete test [5] [16] would be to verify if the 
shear spectral function has a transport peak at to = or not; our calculation addresses this 
question. 

There are many candidate gravity duals of SU(iV c ) Yang-Mills theory. The dual we want 
to apply is improved holographic QCD (IHQCD) as developed in \lT\ [THJ, EH This is a 
bottom-up Anti de Sitter (AdS) gravity +dilaton model which combines elegantly asymptotic 
freedom in the ultraviolet (UV) with a given beta function and confinement in the infrared 
(IR). IHQCD has to be solved numerically and we can so far only compute the imaginary 
part of the Green's function the case k = as a function of to. Since we want to treat the 
entire to, k plane, and for increased transparency, we shall in this note also use a simplified 
version of IHQCD for which the gravity background is analytically known [TS]. It contains 
qualitatively the relevant features of IHQCD and computations in the overlapping region 
confirm that it gives good idea of what one can expect of scale invariance breaking effects on 
Green's functions in SU(iV" c ) Yang-Mills theory. 

IHQCD contains the background metric and the scalar dynamically coupled, and fluctua- 
tion analysis should be carried out taking this into account |21 t l22 | l23]. The transverse tensor 
mode (T\2 with k = (0, 0, k); also called shear mode or scalar mode) does not couple to other 
ones and we thus consider only this case and study the retarded Green's function 

G(to,k;T)= I™ dtd 3 xe wt - lkx3 (ie(t)[T 12 (t,x),T 12 (0,0)]) T = f°° — p(a/ ' fc) , . (l) 

J-oo J-oo 7T to' - tO - l€ 

On scales and notation: in IHQCD, the QCD-like scale A arises as an integration constant 
of the dynamical dilaton field equation, giving the z dependence of the dilaton field, z being 
the fifth dimension. In the model studied here, it is introduced as a constant of dimension 1 



by the z dependence of the ansatz (15). The spin 2 glueball spectrum (squared) then will be 
4A 2 (2 + n), n = 0, 1, 2, .. (analytical), and the transition temperature T c = 0.400A (numerical 



1 A computation in a non-conformal theory but without a phase transition has been carried out in [13] 
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result). Energy (and wave number) can be made dimensionless by scaling with either T or A 
and we shall define 

u = ttTu = 2Aw w 5.00T c w, (2) 
similarly for k. As the dimensionless fifth coordinate we shall use 

y = AV. (3) 

Unless specifically noted, results for the retarded Green's function G are given without writing 
explicitly the dimensionless factor C 3 / (AirG^), where C is the AdS radius and G5 the 5dim 
gravitational constant. 



2 Gravity dual 
2.1 IHQCD 

IHQCD is based on the gravity + scalar action, in standard notation, 
1 



S 



Jd 5 x^g[R-l (<V) 2 + V(<f>)] - J d< 



16vrG 5 

The background, with boundary at z = 0, is chosen to be 



x v — 7 



ds 2 = b 2 (z) 



-f(z)dt 2 + dx 2 + 



dz* 



/(*) 



(4) 



(5) 



where b(z) , f (z) , 4>{z) are solutions of the first order system 



W 
b 
A 
/ 



4bW 2 - 
-b 2 W, 



Uwf + Hv), w = -b/b 2 , 



§ XV bW, 
3/6W, 



(6) 

(7) 
(8) 

(see Appendix A of [E]). The coupling of the boundary theory is associated with X(z) = e^ z \ 
and a specific combination of the fields, arising in connection of solving Einstein's equations 
(b = db/dz), 

/3(A) = ~ (9) 
b/b 

is associated with the beta function of the boundary theory. An important general property of 
the solutions is that b and / are monotonically decreasing and A is monotonically increasing 
as functions of z. 

The solutions have to satisfy the following properties. In the UV, z — > 0, 



b(z) 



C 

z 



(10) 
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Imposing asymptotic freedom to leading order, one then has, 



/3(A) = A = -6 A 2 = -zf-, (11) 
b/b dz 

implying that the leading term of \(z) near the boundary is 

b X(z) = j, L = log~ (12) 
L Az 

A = constant of integration. This is how the dimensional transmutation of SU(A r c ) YM theory 
appears in the gravity dual. Note that on the boundary A — > but </> ~^ ~ 00 • 

Including more powers of A in the beta function, one can derive expansions in A for all the 
fields. For the record, we have collected some of them in Appendix A. 

In the IR, for z — > oo, the crucial confinement criterion is that out of the many different 
solutions of Einstein's equations (Tsl) one should choose the one satisfying 



/3(A) = -l A 



3(0-1) 1 / 1 



4a log A \log 2 A 



(13) 



The parameter a > 1 here is connected with the IR behavior of all the other background 
functions so that, for example, 

b(z) — > e - ^** (z*z) p , 5 = number x A, p = number. (14) 

The appropriate IR relations are summarised in Appendix B. We shall next choose p = — 1 
to extend the form to all z and a = 2 to have a glueball spectrum with m 2 ~ integer }25J . 

2.2 Model: approximate version of IHQCD 

Our model for simulating the essential properties of IHQCD is defined by a simple ansatz 
[E] for b(z): 

C C A 

b(z) = -exp(4 AV) = — exp(-|y), (15) 

where we often set C = 1 and where we use the dimensionless distance variable y in Q. 
Starting from (15) one can derive all the other bulk fields. First, integrating ([7]) leads to 



which integrates to 

A(y)/A = exp(^(y+§) (^ +s Jy~^f\'\ (17) 
Now that X(y) is known, the /3 function as specified by IHQCD is obtained: 



V(y) _ 3 yi + 9/(2y) 
m ~lmJW)--2 X 1 + 3/(2!,) • (18) 
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In the IR this goes like 
/3(A) 



3 X ( , 3 



A 1 + 



;logA 



(19) 



and thus satisfies explicitly the IHQCD confinement criterion. In the UV, A(y) goes to A(0)(1+ 
3y / y/2 + 9y/4 + ...), i.e., this ansatz does not produce a A(y) vanishing logarithmically at the 
boundary. The far UV region is inessential for the treatment of thermodynamics in IHQCD 
|24j and the absence of the " spurious logarithms" |26j simplifies the AdS analysis significantly. 
To get f(z), we note that Q integrates to (C is a constant) 

C/b\z). 



(20) 



The greatest utility of the ansatz (15) lies in that this is integrable in closed form. Using 
/(0) = 1, f(y h ) 

(y - i) e » + 1 y h -i-(y 



= 0, one has 

m = i 



l) e V-Vh 



l)eVh + 1 



y h - 1 + e-Vh 



1 



y h - 1 + e~v* 



n=l 



y h + n - 1 



n< 



(-l)^(y h -yr. 



(21) 



Finally, the scalar potential V(X) can be solved from ^ with the answer (keeping C here to 
show where the dimensions of V come from) 
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£2 



* y [(\y 2 + ly + l)f(y)-(\ +\y)yf(y)]. 



V(\(y)) 

In the vacuum / = 1, /' = and, using (17), the potential at large A is seen to behave as 

V(X(y)) = ^A 4 / 3 v^A(l + 0(1/ log A)), 



(22) 



(23) 



again in agreement with the confining potential criteria in IHQCD. A conceptual problem 



with V(A) in (22) is its dependence on y^. However, what matters for solving the background 
is the range < y < yh, and in this range the potentials computed from (22) for / = 1 and 



f = f(y) are very close to each other. They differ maximally at y = y^ and this difference is 
maximally 4.7% if yh = 0.79. 

2.2.1 Thermodynamics 



To assess the relevance of the model b(z) in (15), we first work out its thermodynamics. 
From AttT = —f(zh) and (21), one obtains the entire thermodynamics parametrised by yh 
(including now the factor C in b(z)): 



T{y h ) 
s(yh) 
p{yh) 



A 



3/2 

Vh 



27r y h - 1 + e-yh 



2 AT, 



4G E 



3/2 



-Vh 



dx [— T (x)]s(x). 



(24) 
(25) 
(26) 



Uh 
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Using (17), we may replace yh by A^, which is used in [18J. The overall structure of T is such 
that it diverges both for small yh (ttT —> A/ \fy~h-, the big black hole limit) and for large yh 
(ttT —t- 2 A-y/yh, the small BH limit). Thus, T(yf l ) has in between a minimum which lies at 

y h = 2.149, T min = 0.3962A. (27) 

For T > T m [ n the system is in a deconfined plasma phase with pressure ~ A^?, but actually 
in a part of this range, T c > T > T m i n , the plasma phase is metastable. The stable phase for 
< T < T c is a thermal glueball gas phase with pressure ~ N®, i.e., in this approach p = 0. 
The transition temperature T c is where the pressure vanishes, which happens at 

yh = y c = 1.6863, T c = 0.4000A = 1.0096T min . (28) 

In the UV, where yh — > and ttT — > A/ \fy~h, 

P £ 3 7T 3 



T 4 4G B 4 



(29) 



independent of the IR scale A. This can be used to fix the dimensionless combination 
C? I (4G5), which then becomes ~ N%. For example, in the present AdS/QCD context one may 
match at some large T 3> T c to the pressure of weakly interacting gluon gas, p/T 4 = tt 2 N^ /45. 
Weakly interacting M = 4 supersymmetric YM matter has p/T 4 = tt 2 N 2 /Q while the strongly 
interacting one has p/T 4 = tt 2 N 2 /8 with 

£ 3 N 2 

(30) 



4vrG 5 2vr 2 ' 

The resulting thermodynamics is plotted in Fig. [T] The interaction measure is broader than 
that of pure YM at large N c and the latent heat smaller. The latent heat is related to the 
large T asymptotics by 



ec_ = sc_ 32 (y c -l + e~ycf „ ycl _ 
T 4 T3 3 y6 6 T 4 



(31) 

r-s>oo 



The smallness of L is due to the fact that the y c -dependent term at y c = 1.6863 is already 
considerably reduced from its maximum value 1/8. There is no parameter with which one 
could tune L to a larger value. 

Finally, the sound velocity is given by 



2 _ 1 3 - y h - (3 + 2y h )e-y* _ 1 1 



^ = 2 ^ ^ 1 ya > = 2 l (-32) 

s 2, + 2y h y h -l + e-Vh 3 ( 3 vrf) 2 

At T c c 2 = 0.024. Note that the Chamblin-Reall background used in [13] has a constant sound 
velocity and does not have a phase transition, since for it always e = constant x p. 
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Figure 1: Thermodynamics from the model b(z). The curves are normalised to their T — > oo asymp- 
totic values, which are ~ N%. At T c there is a weak first order transition with e c /T* — 0.06. The right 
panel shows the region around T c with the supercooled (dashed) and unstable (dotted) phases. The 
curve for pressure is multiplied by 10 for clarity. 



2.2.2 Glueball masses 

Glueball masses are obtained as poles of Green's function and we shall later see how this goes 
through. Effectively, this leads to a Schrodinger-type equation, 

- <(z) + V Schl {z)^(z) = m 2 ij(z), (33) 

where, for the background ^ and for scalar glueballs, 

v - 3 CK P^ + x+J* x- m - (U) 

Vsch <-2{b + 2v) + X +3 bX> X = ^T' (34) 

for tensor glueballs the X— terms are missing. With the above explicit expressions it is trivial 
to evaluate the potentials. Scaling out A by writing V/A 2 — > V, Az — > z, the tensor potential 
is 

15 

V t (z) = ^ + 2 + z 2 , (35) 

and the scalar one 

1 / 2187 \ 

Vs{z) = Vt{z) + (3 + 2*2)2(9 + 2*2)2 (--T - 5346 - 1224 ^ + 21fe4 + 4fe6 J • (M) 

The tensor equation is analytically solvable |25J and leads to the mass spectrum 

ml/ A 2 = 4(n + 2), n = 0,1,2,... (37) 



6 



Solving numerically the Schrodinger equation, one finds the scalar glueball mass spectrum 

m|/A 2 = 6.27, 10.27, 14.25, ... (38) 

The lowest tensor and scalar states are in the ratio niTi/msi = 1-13 (in the full IHQCD this 
ratio is 1.36 [20 1) , but the m 2 spacing is about 4A 2 also for scalars. In the full IHQCD the 
ratio ms2/msi = 1-61 while here it is 1.28. 



3 Fluctuation equation 



To calculate correlators [U El |3l 2] of the operator Tyi in the background ^ one has to solve 
the fluctuation equation 



i + |lo 6 (6V)U(^ 

Here h = hx(z), K = (u, k), and 

d 4 K 



h + Ph + Qh = 0. 



h(x, z) 



(2tt)< 



h (K)h K (z), h K (0) 



(39) 



(40) 



is the fluctuation of the 12 component of the background metric: gi2 = b 2 (z)(l + h(t, x 3 , z)). 
This tensor fluctuation does not mix with other ones. In IHQCD, the scalar one does and the 
fluctuation equation permitting computations of the Green's function of the scalar operator 
F 2 „ will contain more terms [2T1 l22l 123] . 

The characteristic exponents of (39) at z = z^ are ±iu)/fh. Eq. (39) should be solved so 
that at the horizon (infalling boundary condition, di are calculable coefficients) 

h K (z->z h ) = (z-z h r/f»[l + d 1 (z-z h ) + d 2 {z-z h ) 2 + ...] (41) 

and that at the boundary, hx(0) = 1 (correct normalisation). These solutions will satisfy 

h- K (z) = h* K (z). (42) 

Given a solution, a computation along the lines of [27] shows that the retarded Green's 
function is obtained from 



G{K) 



1 



16vrG 5 



fb s hxh-K + (counter terms) x hxh-K 



(43) 



where the hxh-x = hxh* K terms with unspecified coefficient are real. 

For the imaginary part of G, the spectral function, the situation is particularly simple: 



ImG(fc) = p(u,k) 



1 

16vrG 5 



fb 3 Im h K h_ K 



1 

16ttG 5 



fb 3 W(h K ,h* K )/(2i) 



(44) 
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where one has noted that hx and h* K are two independent solutions of (39) and that their 
Wronskian arises when ImG is evaluated. However, the Wronskian of any two linearly inde- 
pendent solutions of (39) is, integrating W/W = —P, 

W (K) 



W((pi, 2 ) = 0102 ~ 0201 



63/ 



(45) 



where Wo is a z independent constant (but will depend on K 
is a ratio of two Wronskians and independent of z. 

To find the constant Wo for the Wronskian W{hx, h* K 
z — > Zh using the expansion (41) and find 

2iu , o . 



oj, k). Thus the result in (J44J) 
we simply evaluate it in the limit 



W(h K ,h% 



K) 



fb 3 



The result thus automatically contains the entropy density, which in IHQCD is 

A 1 



s(T) 



-b 6 {z h ). 



(46) 



(47) 



4G 5 4G 5 

However, so far this is for a solution normalised as in ( |41[ ), which leads to some value /ik(0) 
at the boundary. Correcting for the normalisation, the result is 



p{ ^ k) = i s{T) \h^w- 



(48) 



A practical note on the application of this formula: in it hx(z) arises from numerical 
integration which cannot be extended to z = since this is a singular point of the equation. 
However, in (48) hjc(z) has to be evaluated precisely at z = and this value has to be found 
by extrapolation. Numerical effect is significant at large u, where |/i#-(0)| ~ 1/oj 3 ^ 2 . 

We shall also need the real part, for example, the static Green's function G(Q, k) is purely 
real. For the real part one has to make further assumptions, corresponding to subtractions in 
the dispersion relation Q. The standard procedure here [3] is to expand the solution hx(z) 
computed with the boundary condition (41) in terms of the two independent solutions at 
z = 0: 

h K (z) = A{K)<f> u (z, K) + B{K)(t> n {z, K), (49) 

where 

4> n = Z 4 (l + C\Z 2 + C 4 Z 4 + ....), 

M = 1 + hz 2 + hz & + ... + c\ogZ(j> n {z,K). (50) 
In IHQCD the expansion coefficients will contain logarithmic UV coefficients of the type in 



Appendix A [T7]. Noting that the normalisation factor ^(0) in (48) is A(K) the result (43) 
can be evaluated using 



hVh 



K K 



1 



^IMO)! 1 



(A*0 u + 5*0„)(.40 u + £0„) 
i + 4clog z + c + 2f 2 + 4||^| + 0{z 2 ) 



z*\A? 

■y2 



uj 2 - k 2 (cj 2 -k 2 ) 2 , A , , B{K) ^,2, 

(41ogz-l) + 4-)-( + 0(z 2 ), 



2z- 



16 



A(K) 



(51) 
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where the last form applies for conformal expansions (50). The scheme is to neglect the real 
divergent terms + constant in (51) and to evaluate the entire Green's function from 



G{K) 



£ 3 B(K) B(K) W{h K ,<j> u ) 



AttG 5 A(KY A(K) W(h K ,4> n ) 



h K 4>n ~ (t>nhx 



(52) 



Equivalently, one may say that this is a way to derive the counter terms in (|43j) [28]. In [5] 
this scheme was tested by evaluating the purely real static Green's function G(0, k) also by 
an ^-integral over the complex G(u,k). It also gave correctly the real and imaginary parts 
on the light cone, for G(k,k). 



Since G(K) in (52) is given as a ratio of two Wronskians, it is independent at what z it is 
evaluated, to the extent hx, 4>u, <Pn are correct solutions. Since hx is integrated numerically, 
it certainly is a solution. In practice the accuracy and constancy then depends on how many 
terms are included in small- z expansions of 4> u ,(j) n . If the numerical integration of hx from 



(39) is extended down to some small e, the ratio of Wronskians can be evaluated at this e. This 
is to be contrasted with (48) in which hx(0) has to be extrapolated to z = for accuracy. 



Of course, (48) only gives the imaginary part while both the real and imaginary parts are 



obtained from (52). 



3.1 Equations for the model 



Next, we specialize to the model defined in subsection 2.2 In terms of y = A z the fluctuation 



equation (39) becomes 



f K { y ) + P{y)Mv) + Q(y)My) = o. 



where 



P(y) = — log(6 3 /^) 



dy 



1 + 



y 



Q(y) 



CO 2 



y-l + il-y^eVh-y' 
k 2 



yf 2 {y) yfivY 

where the dimensionless Cj, k were defined in &2h. 



(53) 

(54) 
(55) 



The equation (53) is to be solved so that at the horizon (the y — yh term is for k = 0) 

+ 



<Px(y) = (y - yh) p 



i+ 2yh{l - 2p) (y-y h ) + 0((y-y h )) 



where 



P 



-IU! UJ 

l- 



4-kT 



yhf'ivh 



Vh 



(Vh - 1 + e 



(56) 
(57) 



For correct normalisation, the solution is divided by 0^(0). 
The Green's function now is 



G(K) 



2A 
16vrG 5 



fb^^/y (f>' K (j)-K + counter terms • §k4>-k 



(58) 



9 



For the spectral function one again needs the Wronskian of the two independent solutions at 
the horizon, 



K) 



2iQo 



b 3 (Vh 



(59) 



and the result (|48|) is unchanged. For the real part and evaluation of the counter terms in 

(60) 



(58) the expansions around y = have to be worked out. One writes 

<j> K {y) = A(K)MV, K) + B(K)4> n (y, K), 



where 



4> n (y,K) = y 2 
1 



l + l(2 + k 2 -u 2 )y + ... 



(61) 



<f> u (y, K) = 1 + {uj 2 - k 2 )y - - [{uj 2 - k 2 ) 2 - (Co 2 - k 2 )} log y ■ Mv) + 0(y 3 ) + (62) 



(note no y 2 term; it is only in (fi n ). The real counter terms neglected now are 



2A 4 



UJ 



1 



y 2 

and, solving A and B with Wronskians 



[ ( ^_^_(^_^)](21ogy+l) 



C 3 

G(K) = ■ A 



A(KY 



4ttG 5 

Now the correct dimension 4 comes from A 4 . In units of ttT, using (24): 



G(K) C? 16(y h - 1 + e-^) 4 B(K) 
(ttT) 4 ~ 47^5 y\ A(K) • 



(63) 



(64) 



(65) 



We are interested in the high T phase and, as was derived earlier, this means T > T c 
0.4000A or y h < y c = 1.6863. 

3.2 Vacuum spectral function and the limit at large u or k 



To obtain the vacuum Green's function of the model we put / = 1 in (53) and obtain the 
equation 

y<t>"{y) + (-1 - y)<f/(y) + K 2 <P(y) = 0, K 2 = G? — k 2 . (66) 



This also gives the leading terms of the uj, k ^> A, T limit of the solutions of ( 53 ). Understand- 



ing this limit analytically is important for numerical accuracy. This equation is the equation 
for the confluent hypergeometric function with a solution 4>i growing like a power at large y: 

Mv) = V 2 U(2 - K 2 , 3, y) -»■ / 2 [1 + K 2 (2 - K 2 )y- 1 + 0(y- 2 )] . (67) 
Its small y expansion is 



i(y) 



r(2 - k 2 ) 



l + k 2 y + -k 2 (l-k 2 )[\ogy + i)(2-k 2 ) + 2 lE -^ )y 2 + 0(y 3 ' 
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The second solution 02 (y) 



,.2r(2) 



(v) 



e^y ^ +1 + .., grows exponentially and is un- 



acceptable. If we now decompose the solution <pi as in (60), B is immediately read as the 
coefficient of the y 2 term and A = 1/T(2 — X 2 ). The vacuum (this is vacuum since no T is 
involved) Green's function then is 



Gvac(K) = A 4 -\k 2 (1 



K 



^(2 - K 2 ) + 2 lE 



Since 



m=0 



1 

m + 1 



1 

771 + Z 



the result in (69) has poles at 



4A 2 (m + 2). 



(69) 



(70) 



(71) 



This reproduces the tensor glueball mass spectrum in (37), now derived as poles of Green's 
function. Taking u — > ui + ie and plotting ImG vac from (69) for finite but small e one finds 
Gaussian peaks at glueball masses. In the limit e — > 



Im G v 



TT 

32 



(u 2 -k 2 )(cj 2 -k 2 



4A 2 



m=0 



771 + 2 



Co'" 



A- 



4A 2 



(72) 



We shall later use this as a baseline to assess the magnitude of finite T effects (see Eq. (79)). 

The above was for vacuum, / = 1. The asymptotic large u, k limit for / ^ 1 can be argued 
as follows. Taking first uj > k > 0, the large K 2 > limit of tp in (69) is 



■0(2 — K 2 ) = 7r cot(7rK 2 ) + log K 2 
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2K 2 12K 4 



+ ... . 



(73) 



Note now that 7rcot(x + iy) —> ^iir if y — > ±oo. If now the analytic continuation is performed 
so that also Im K 2 — > oo, then (69) and (73) give 



ImG(K) _ it 
A4 ~ 2 



^ -k 2 ) 2 -k 2 ) + 0{\) 



(74) 



Both terms are very clearly seen in numerics, but we have been unable to fit a constant (or 
possibly a logarithmic) term nor compute it analytically. 

In contrast, in the static uj = case, K 2 = —k , the constant term (and higher corrections 
0(l/k 2 )) can be evaluated. Now Eq. (69) is purely real. Expanding -0(2 + k 2 ) — > log A; 2 + 
3/(2fc 2 ) - 13/(12/c 4 ) + C(l/A; 6 ) in ((69]) gives a constant term -5/24. However, this is not the 
entire constant term, which is crucial for numerical analysis. A careful analysis (see Appendix 
[C] ) permits one to compute it with the result (116), which is the same as writing (69) in the 
form 



B 
A 



1 



A: 2 (l + P 



V>(2 + A; 2 ) + 2 1E 



-Vh 



+ 



20(y ft 



1 + e~Vh] 



(75) 



Thus there is an overall constant term —5/24 + e y,1 /(20(y/ l — 1 + e Vh )) 
This is clearly seen in numerics. 



-5/24 + T c 4 /(5T) 4 . 
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Figure 2: The spectral function lax G (us, 0; T)/(ttT) 4 in units of £ 3 /(4ttG 5 ) plotted vs w/(ttT) (left 



panel) or divided by uj/(t:T) (right panel), for the model (15). The dashed line is the leading large-w 
limit. 
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panel) or divided by lo/(ttT) (right panel) in exact numerics of IHQCD. 



In the thermal theory with / ^ 1 the glueballs form a glueball gas with pressure 0(N®) 
and AdS/QCD duality says little about its properties pD], it is the phase at p = in Fig. [TJ 
In the high T phase at T > T m [ n the delta function peaks are not seen but still some structure 
remains (see Figs. [7] and [8] below). The broadening of hadronic states in a thermal ensemble 
has been studied extensively in connection with quarkonium physics 



4 Numerical results 

We are interested in computing the finite temperature retarded Green's function of T\i in a 
QCD-like strongly interacting theory, for various frequences and wave numbers k = (0, 0, k). 
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Figure 4: The imaginary time Green's function G(r,k — 0) computed from (77 1. The right panel 
shows the region near r = 1/(2T) on a linear scale, the dashed curves come from full numerical 
evaluation of IHQCD. 



At large T, say, T > 10T C one expects the theory to approach a conformal theory with 
Green's functions explicitly plotted in [U El [8]. In the range T c to 10T C there are important 
nonconformal effects, reflected in the large value of the interaction measure (e — 3p)/T 4 . This 
is the range of T to be studied here. 

We remind that the results for G or p = ImG are given without an overall factor of 
C 3 /(AttG 5 ). 



4.1 Spectral function for k = 

This is clearly a very interesting range due to its relation to viscosity and lattice determina- 



tions thereof [331 [Ml E51 [361 E7J [39] . Consider first the results for the model (15). As discussed 
in subsection |3.2[ the large-w limit is, measuring both G and u in units of ttT, 



Pas(u},k) ImG as (w, k;T) 



7T 

32 



(u> 2 



k 



2\2 



25.02* 2 



k 2 



6(oj - k). (76) 



This is an asymptotic formula and both terms are clearly seen in numerics. However, we 
do not know the terms of lower order in oj 2 — k 2 , but (76) should be applicable only for 

2 > k 2 



25.0T 2 , it is negative at smaller u. 



Results for the model ( 15 ) are shown in Fig.[2]together with the leading asymptotic behavior 



TTUJ 



/32. The spectral function is shown both as such or divided by uj so that the value at 
uj = is essentially the shear viscosity. One observes a sizeable effect near T c but already at 
2T C the curves are very close to the conformal situation, approximated by the 10T C curve. 



The main T-dependent effect is that contained analytically in the asymptotic behavior (76). 
Even for T ~ T c the nonconformal effects disappear at uj>3T c ; remember that the lightest 
scalar glueball mass is ~ 6.2T C . 
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Figure 5: The weighted spectral distribution, integrand of (77), at the middle point r = 1/(2T) for 



both the model (15) (left) and full numerics of IHQCD (right), at various temperatures. 



Consider then the results, shown in Fig. \m computed from full numeric g^of IHQCD, using 
the formula (48). This gives only the imaginary part; due to the complicated logarithmic 



structure of IHQCD (see Appendix A) it will be more complicated to evaluate the Wronskians 
in (52) and compute the real part also. The overall pattern at k = is very similar. An 



interesting detail is that in IHQCD the approach to the asymptotic limit is rather slow, 
possibly related to the UV logarithms. 

The imaginary time Tyi correlator at k = can now be integrated from 



G(r,0;T) 



du cosh[(l-2Tr)f a;] 
— p(u, 0; T) . ? 



(77) 



with everything in units of ttT. Note first that in the J\f = 4 conformal case the answer is |8j 



G(r,0;T) 



47T 5 



[a5,TT) + ((5,l-Tr)} + Gp(TT), 



(78) 



where the first terms come from the integration of 7T£j 4 /32 and Gp from p(ui,0) — 7rw 4 /32. 
The thermal part Gp varies very slowly from Gg(| ) = 0.0710 at the symmetry point to 
the value 0.13 at r = 0. At TV = \ the vacuum part, i.e., the £ function terms, contribute 
93£(5)/(27r 5 ) = 0.1576 and dominate totally at smaller t. 

The outcome in the nonconformal cases is shown in Fig. [4j The T = 10T C curve is very 
close to the conformal curve (78). A particularly interesting point is the middle one, r = 



1/(2T). Fig. [5] shows the corresponding integrand of (77), in which large w values are strongly 
suppressed, Fig. [6] shows the resulting value of G(t = 1/(2T)). 

To assess the significance of the result, it is very useful to generate a comparison curve 



by inserting the vacuum spectral function (72) with delta function peaks at the position of 



2 The parameters of the scalar potential were otherwise the same as those in [20] but Vb 
Vo = 0.04128; this makes numerics simpler and does not affect the results. 



1 instead of 
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Figure 6: The temperature dependence of G(1/(2T)) for both the model (15) and full numerics of 
IHQCD. The arrow shows the value for the conformal j\f = 4 theory [8; The dashed curve shows 
the plot of (79 1, generated from the T independent vacuum spectral function ( 72 1 via the spectral 
representation (77). 



glueball masses to the representation (77) at r = 1/(2T). The answer is (A = 2.5T C 



G vac (r = 1/(2T)) _ 1 {2.5T C \ 5 v^T+2(m + l) 



(7rT) 5 2\itTJ ^ Q s■mll{2.5^/m + 2T c /T)' 

The resulting curve is also plotted in Fig. § the limit at large T is 93C(5)/(2tt 5 ), the same 



as for the vacuum curve in (78). One sees that the imaginary time correlator generated by 
the vacuum spectral function behaves similarly to the genuine finite T correlator, finite T 
just enhances the magnitude. The vacuum correlator has nothing like a transport peak at 
uj = 0, it is identically zero up to the lowest glueball mass at 2\/2A. Still the outcomes are 
very similar. 

We do not attempt a detailed comparison with lattice Monte Carlo data. However, the 
results are in qualitative agreement with the plots in [35]. An agreement confirmed by careful 
analysis would support the present picture of strongly interacting QCD matter, but Fig. [6] 
illustrates the difficulty of making definite conclusions. 

Since one knows that the viscosity prediction holds independent of T it is useful to present 
the data so that this is explicitly seen. Data for the dimensionless quantity 

(80) 

s(T) ui 

computed for the model (15) in the region of small uj, are presented in Fig. [7J left panel. The 
normalisation is such that the standard viscosity prediction r]/s = l/(4.7r) corresponds to the 
value 1 at uj = 0. The right panel of Fig. [7] shows the spectral function with the asymptotic 



behavior (76) subtracted. The subtraction is reliable only for uj>5.0T c /(ttT). 

Finally, Fig. [8] shows the subtracted Green's function plotted in non-thermal units, vs. 



uj/(2A), computed for the model (15). This is what remains of the infinite set of delta function 
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4 7T p(u, 0; T) - p a 




Figure 7: The imaginary part of the Green's function G(w, 0;T) normalised as in (80 1 plotted vs 
w / (ttT 1 ) for small w (left panel) or larger u) with the asymptotic behavior ( 76 1 subtracted (right 
panel). The subtraction is reliable only for uj/(ttT)>5.0T c /(ttT). The model (15) is used here. 




4tt p(co, 0;7)-p a 




V U>/(2A) 

5T 8 ^-40 12 



Figure 8: The imaginary part of the Green's function G(w, 0; T) with the asymptotic behavior sub- 
tracted, normalised as in (80) plotted vs Co — cj/(2A). The model (15) is used here. 



peaks in (72). One may try to define a glueball mass at finite temperature by the position of 



the first peak in this curve. This grows very accurately linearly with T: 

m G (T) = — ^ = 1.016- - 0.241, 
mG/v2 J-c 



(81) 



where the lowest tensor glueball mass is vtlq = 



16 




4.2 Real static correlator, oj = 

The static correlator is defined by 

r/3 r ■ 

'0 JO *' H ~ r Jo 

and can also be related to the w-dependent spectral function by 



G(r)= dT(Ti 2 (T,x)Ti a (0,O))r= / dr G(r, x) = — =- / dk ksm(rk)G(k) 
Jo Jo r Jo 



G(k) 



du 2p(co, k) 



7T 



(82) 



(83) 



It is determined numerically by solving (39) with co = so that the boundary condition at 



the horizon is (41) with oj = and solving G(0,k) from (52) as a ratio of two Wronskians 



It is also important to eliminate the large- A; behavior given in (75), including the constant 



term. Upon transforming to coordinate space these give rise to 5(x) or its derivatives. To get 
G(0, k; T)/T A one further has to use (65). 



Subtracting the expansion (75) from the numerically computed B/A gives the result in 
Fig. [9] for the static correlator. The curves are of the same shape as those in the conformal 
case jS], where one plots vs k/(7rT). When T — > oo the value at k = of the curves approach 
the value — 7r 4 /10 of the conformal case. The dominant scale breaking effect is the rapid 
decrease of the correlator when T decreases towards T c . The next subsection explains how 



this is via Eq. (83) related to the behavior of p(u, k) near the light cone, see Fig. 12 



In the conformal case one can define a second order dissipative coefficient by the expansion 



G(0, k; T) = G + \ n{T)k 2 + 0(k 4 



and its value is 140 



£ 3 



AttG 



- 3 7T T 



I nJ.t 2 



4vr 2 T' 



(84) 



(85) 
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Figure 10: The expansion coefficient k(T) defined by (10 1, normalised to the value of k(T)/T 2 at 

iot;. 



Inserting the normalisation (30), our computation of k at 10T C is 1% from the conformal 
value in (85). At smaller T there is additional suppression on top of the T 2 dependence. This 



is quantitatively shown in Fig. 10 



Given the momentum dependent static Green's function G(k), the coordinate dependent 
correlator G(r) can be computed by carrying out the integral p2j ). The shape of the curves 
is very similar to that plotted explicitly for the conformal case in [8]. 



4.3 Green's function for any u, k 

The dominant feature on the u, k plane is the structure around the light cone. In [8j this was 
analysed analytically and, for example, one could prove that exactly on the light cone 



G(k,k) = (l + iV3)- 



r(i/3) 



6V3r(2/3) 



&4/3 + Q ^ ra (0.136092 + i 0.235718) fc 4/3 , k = k/(*T). 

(86) 



Results for the spectral function as a function of uj at fixed values of k are shown in Fig. [TT| 
The left panel shows what effects changing T at fixed value of k/ijrT) has: the correlation 
function monotonically decreases. The right panel shows that essentially the conformal set of 
curves is obtained for T = 10T C , for the model (15). 



Another useful way of analysing the light cone region is to see how the spectral function 
varies when crossing the light cone perpendicular to it, i.e. at fixed <jJ + = (oj + k)/\/2. The 



outcome is shown in Fig. 12 at a fixed value of uj + . At larger uj + the peak gets narrower 
~ 1/Zc 1 / 3 but higher ~ /c 1//3 so that the integral over oj, which via Eq. (83) leads to the static 



correlator, is essentially constant. This is as shown in [8] for the conformal case. The new 
element is that when T decreases towards T c , the correlator also decreases. After integration 



over uj according to (83 ) this correlates with the rapid decrease of the static correlator, already 
shown in Fig. |9] 
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Figure 11: The imaginary part of the Green's function G(uj, k = 2-kT; T)/T 4 in units of £ 3 /(47rG 5 ) 
plotted vs lj/(itT), for the model (15 1. The right panel shows the behavior near the peaks at k = 
2, 5, IOttT. The curve 0.236a; 1 / 3 is the conformal behavior of the correlator on the light cone, U) = k. 




Figure 12: The quantity (p(u), k; T) — p as (tv, k))/ui when crossing the light cone in the perpendicular 
direction, i.e., at fixed w+ = {u-\-k)/\2 (here uj + — 10-\/2) with varying for various temperatures. 

5 Conclusions 

We have in this paper computed the retarded Green's function G(ui, k; T), k = (0, 0, k), of the 
1, 2 component of the energy momentum tensor in the high temperature phase of a QCD-likc 
strongly coupled theory with a first order phase transition and glueballs, using gauge/gravity 
duality. Temperatures between T c and 10T C were numerically studied. Results were given for 
the uj dependence of the spectral function at k = 0, from which one could compute the T 
dependence of the imaginary time correlator G(r; T), measured on the lattice. The real static 
correlator at u = was also computed and related to structure observed near the light cone, 
uj ~ k. Overall, at 10T C the correlators were very close to the conformal situation [8], but 
large effects were observed when T decreased towards T c . 

The gravity background used here is IHQCD as developed in [E]-[20] for k = and a 
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simplified and quasianalytically tractable version thereof for all u, k. IHQCD describes very 
successfully the properties of large N c Yang-Mills theory both at T = and at finite T, but its 
gravity background, in general, has to be evaluated numerically. The model used here retains 
the attractive features of IHQCD in the infrared, so that it has confinement, glueballs and a 
first order phase transition, while simplifying the complicated logarithmic structure of IHQCD 
in the ultraviolet. The chief virtue of the model is that both b(z) and f(z), appearing in the 
fluctuation equation, are analytically known. Another analytic approximation to IHQCD 
\2A\ 141] , would be quantitatively more accurate but does not permit analytic evaluation of 
f(z). In any case, we expect that the results of this paper should be qualitatively valid for 
IHQCD, too, for all uj, k. This is confirmed by numerical computations at k = 0. 

There is ample phenomenological evidence for the strongly interacting nature of QCD 
matter at T>T C (p = |pideab small viscosity, jet quenching), but one has not been able to 
develop a first-principle numerically verifiable criterion for this. We have suggested that a 
quantitative verification of our predictions could serve this purpose. 

There are several obvious directions of further study. For the first, a full numerical compu- 
tation in IHQCD should be carried out, also for nonzero k. This is technically complicated 
due to the fact that the expansions near the boundary are not power series in z but contain 
powers of log z, studied here in Appendix A. Secondly, to extend the calculation to correlators 
of the scalar operator Tr F 2 U one should use the more complicated fluctuation equation in 
[23], coupling the dilaton and the scalar component of the tensor fluctuation. A third open 
problem is the computation of more terms in the large u expansion (74). Finally, it will be 
interesting to compare the results with those from next-to- leading order perturbative QCD, 
when they are available. 
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A Appendix: UV expansion in IHQCD 

In this appendix, we shall give small- z logarithmic expansions of b(z), X(z) and W(z) = —b/b 2 , 
adding a few terms to those given in Appendix A of [T7j. Also f(z) will obtain similar 
logarithmic corrections to the z 4 term by integration of (20). 
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The expansions are given as powers of X(z) or as inverse powers of 

1 



L = lo£ 



Az'' 



(87) 



where A is the integration constant of 



dX 

dz 



-fiWb(z). 



The correction terms arise from the UV beta function expansion written in the form 



b P(X) = -(6 A) 2 -^(M) 3 -3(6oA) 4 -S(M) 5 + .., (89) 

°0 °0 °0 

where only 6q, b\ are scheme independent. It will appear that all quantities are functions of 



b2, 



6nA, 



6 i+l' 



6 2 ' 

°0 



(90) 



(b = bi/bQ is a standard notation, not to be confused with b(z)). 

First, by integrating db/b = dX/(3(X) = dboX/ \bQ(3(X)), (note the different constants &o an d 
bo) 



b 



cxp 



1 

k>X 



(6oA) 



bi/b 2 



1 + 



'>■! ,2 



°0 



6 A + 



6(l + ^-2^) + | + ^ 

°0 °0 °0 



1 



The constant of integration bo is fixed so that 6oA appears in the powerlike term. 
For W = —b/b 2 one can derive 



W(X) = W(0) exp 



A 2 J 



dX 1 ^- ) , W(0) 



1 

£' 



(b X) 2 + . 
(91) 



(92) 



Expanding this one has 



CW = 1 + - 6 A + 



8 2 6 
81 



The z dependence is obtained by combining the equations W = —b/b 2 and (3 = bX/b into 

dX 



(88) so that, integrating, 



dz 



-pWb' 



(94) 



Noting that W and 6 contain the constants 1/C and 6o this equation naturally defines the 
scale 

bo 



A 



(95) 
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of z. Eq. (94) now has to be inverted for A: one expands the RHS around A = using the 
expansions above, integrates term by term and takes the limit A — > 0. Finally, the result is 
inverted recursively to obtain A = \(z). The expansion one obtains is 



6nA 



1 61ogL 



L 



L 2 



+ & 2 log 2 L-6 2 logL + 



a M 1 



+ 
+ 

+■ 



-b 3 log 3 L + ^b 3 log 2 L 



3b 



b 3 o 



2b 3 + -b logL 



4 2 , 1 ,o 6 3 
9 + 3 & "2 6 + 2l| 



L 4 



+ 



6 4 log 4 L 



13 



S 



b" log' 5 L + - 6 + - 6 + 66" 



104 
81~ 



9 6 27 6g 



,62 
5^1 



log 2 L + 46" 



. . , 1 , b-x 1 64 
t \h i — A - h— -I - 
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logL 
(96) 



To get corrections to b(z) one writes (88) in the form 



IdA 

2 dL ' 



(97) 



uses (96) and expands in A. This leads to 

b{z) 



-4 



C 

z 
C 

z 

710 



4 

1 - n M + 
y 



44 2 l| /l ^2 

6 (60A) 

81 9 ' 



/2408 80 L 4 6 2 \ , , , ., , 
(2187-81 6+ 27 63) (M) + ° (A 



4 2(1861ogL + 22-96) 
1 ~9L + 81L2 + 



(98) 



20, 



2187 81 



b 2 + 



4 bo 



+ Trrb 



27 b 3 Q V 81 9 



logL + ^& 2 log 2 L 



The three first terms in the second form were given in (A. 28) of [17] . 

Finally, the UV expansion of the potential V is very simply directly obtained by inserting 
the expansions of W and /3 to V(X) = 12PF 2 (A)[1 - (/3/(3A)) 2 ]: 

'23 4 



V 



12 

£2 



1 + - b X + 
y 



(99) 



B Appendix: IR expansion 

One cannot uniquely solve the equations Q at large z, large A, but assume that b(z) behaves 
as 

b( z ) = b e-^ a (Ez) p , a>l, p real. (100) 

E here in the IR is a number times A in the UV and, to emphasize this difference, we 
use a different letter for it. Inserting this to the equations of motion above and evaluating 
successively b/b,bW,W,W /W one can write <j> in a form which is easily integrable. If Aq is 
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the constant of integration of this equation, one finds that, up to corrections 1 + 0(l/z a ) 
l + 0(l/logA), 



b(z) 

W(z) 
X(z) 

/3(A) 
V(z) 



b e-^ a > 



Z>Z) 



a JEz)° 



^(llogA)(^)/(-), 
)«-i-p= i A 2 / 3 a(| logA)^- 1 )/^), 



60 

A ' 



3(Ez) a /2 



(s,) 3 ( 



a-l-2p)/4 



(B*)- = A 2/3 (Az)P - a /2+l/2 5 ^ = ( 2 ]ogA )l/« ) 



■ A 1 + 



ft 



9W 2 (z) 



2a(Ez) a 

?^ e 2(H«)«^ s ^2(Q-l-/-) 



lAfl + ^W ' 



4a log A log 2 A ' 



^A^flogA)^- 1 )/^), 



(101) 
(102) 

(103) 
, (104) 
(105) 



where always A = A/Ao, apart from the leading term in /3(A), where the leading term is — 3A/2 
for any Ao- By choosing p suitably one can make any of the quantities a pure exponential 
in z a . However, when expressed in terms of A, the leading behavior is always independent 
of p. Note that 6A 2 / 3 = b s ~ (log A)^ -1 ^ 2 "* 1 is the string frame b, which is used to express 
the confinement criterion. Note also that the confining beta function is not assumed here, it 
follows from the assumed form of b(z). The scalar potential V satisfies V = 12W 2 (1 — X 2 ), 
X = /3/(3A). 



Finally, Eq. (20) now is 



f(z) = Ce 3 ^ a (EzY 



3p 



(106) 



In our model with a = 2, p = — 1 this is applies at all z and is simply integrable in closed 
form. If one wants just a large z approximation for f{z) one can integrate it up to corrections 
of order l/z a by extending the integration to z = and imposing the condition /(0) = 1. 
This leads to 



/(*) 



1— a— 3p 



= 3S Q ( Z -- 2 «) 



4irT = 3aH (Hz/i) 



a-l 



(107) 



C Ultraviolet limit of the static correlator 

In this appendix we derive the high-fc limit the static correlator G(u = 0,k;T). Eliminating 



the first derivative in ( 53 ) by writing 



(108) 



and setting uj = 0, we have the equation 
k 2 ( f'(v)\ 2 1 i 



i>"{y) 



(f(y)\ 
yf(y) \V{y)J 



+ 



4 2y 4y 



iP(y) = [k 2 F(y) + G(y)] ^(y) = (109) 
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We wish to solve this around y = for k — >• oo. The systematic approach given in [30J starts 
by defining 

* - f -7^=, = f ^ V2 (no) 

7o ^xf{x) \dy J 

The idea with replacing y by the new variable £ is that the resulting equation is in a 
Schr6dinger-like form in which only k 2 multiplies the function to be solved. Introducing 
W eliminates the first derivative generated. The solution then is given in the forr^ 



where the coefficients A s and B s are computed from the corresponding expansion coefficients 
A s describing the solution away from the singular point £ = 0: 

2s 2s+l 

mo = J2 a ^~ jA **-i(0, b.(s) = «i(2)r i-1 A2 S -i-i(e), (112) 

j=0 j=0 

where aj{v) = fR*^!/ 2 — ~~ l) 2 )/8fc. Functions A s in turn are given by 
A s+ i(0 = -~A' S (£) + \J d^(0M0 

2 2 J Vyfiv) 

The integration constants in each step are chosen so that when we retrace our steps back to 
the original mode function <j>(y), we have 0(0) finite and equal to 1. The function ^(y) in 
the above recursion relation is given by the coefficient functions F(y), G(y) of eq. ( 109| ), and 
reads 

Working out the series expansions of A s , s = . . . 4 (with A s =constant), we can read out 
the coefficients of the y° and y 2 terms in 

Hy) = (b 3 fVyr 1/2 ^(y) = i^VIr^MO (us) 

to finally arrive at 

B W t ■ 3\ r o Ac X 5 



fc lnfc + 7_E r In k + 1e) + — ; r • (116) 

A V 4 7 \ ' J 24 20(y h - 1 + e-v») v ' 



3 For the coefficients A s and B s , we follow here Exercise 5.2 of Chapter 12 in [30] 
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